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1. INTRODUCTION 

Today's power industry needs the development of more complex nonlinear power system models 
and optimization techniques to solve them, these are called the optimal power flow problem (OPF) 
techniques. OPF is one of the most important tools for inefficient planning and controlling the operation of 
power systems. It was first introduced by [1]. The OPF procedure consists in choosing the optimal values of 
the control variables of an electrical system to optimize an objective function while satisfying the constraints 
of equality and inequality of the system [2]. Several objective functions related to the electrical system can be 
optimized, such as: minimize total generation cost (fuel cost, wind energy, cost of flexible transmission 
system (FACTS) cost, etc.), transmission losses, voltage deviation, voltage stability index, toxic gas 
emission, system safety, etc. [3-5]. The OPF problem can be considered as a large problem of nonlinear 
optimization with constraints. The optimization problem solved by several developed mathematical 
techniques, these techniques may be classified into two groups; conventional methods and recent intelligence 
methods (evolutionary or metaheuristic methods). Recently, several evolutionary or metaheuristic 
optimization methods have been proposed to get the best solution to the OPF problem. 

Metaheuristic algorithms (MAs) mark a great revolution in the field of optimization, allow finding 
one or more solutions to complex optimization problems [6]. According to [7], the MAs can be regrouped 
into four main categories: evolution-based methods, physics-based methods, human-based methods, and 
swarm-based methods. Several metaheuristic algorithms are implemented in electrical power system for 
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solving the optimal power flow problem with different objective functions such as moth flame optimizer 
(MFO) [4], enriched brain storm optimization (EBSO) [8], moth swarm algorithm (MSA) [9], particle swarm 
optimization (PSO) [10], cat swarm optimization (CSO) [11], chaotic whale optimization algorithm (IABC) 
[12], improved strength Pareto evolutionary algorithm (IEA) [13], symbiotic organisms search algorithm 
(SOSA) [14], stud krill herd algorithm (SKH) [15], modified Grey wolf optimizer (MGWO) [16], differential 
search algorithm (DSA) [17] and integrated algorithm (IA) [18]. 

This paper presents one of the newest flexibility and efficient optimization metaheuristic method, 
called ant lion optimization (ALO). Recently, many researchers are interested in this method for solving the 
optimization problem, as in [19-20]. In this study, the proposed method has been applied for solving the OPF 
problem for large scale power systems which is the Algerian 114-bus power system. Two different cases are 
considered, with and without the presence of SVC devices. The objective function used in this paper is 
minimizing the total fuel cost (TFC). 


2. MODELING OF SVC DEVICE 
The static VAR compensator SVC is modeled by shunt variable admittance. Since the power loss of 
the SVC device is assumed negligible, so the admittance is assumed purely imaginary as follow: 


Ysve = Jbsye (1) 


The susceptance boy- can be capacitive or inductive to respectively provide or absorb the reactive 
power Qcsyc. The placement of SVC devices in this study is installed in the power system as a PV bus with 
the real power generation equal to 0 MW. The reactive power Qsyc absorbed by the SVC device and also 
injected into node i is given by (2): 


Qsvc = WV? Deve (2) 


3. OPTIMAL POWER FLOW (OPF) PROBLEM FORMULATION 
3.1. Formulation problem 

The solution of the OPF problem aims to minimize or maximize an objective function for getting an 
optimal adjustment of control variables in the power system by satisfying both constraints, equality and 
inequality constraints. Generally, the optimization problem can be represented mathematically as follows: 


Min F(x,u) (3) 
Subjected to g(x,u)=0 (4) 
h(x, u)<0 (5) 


where: F represents the objective function, x represents the vector of the state variables and u represents the 
vector of the control variables. 


3.2. Objective function 
The objective function in this study is the quadratic equation of generation fuel cost of each 
available conventional generator subject to operating constraints and formulated as follows: 


NG 
C, ( P5:) = D4, + Po +.6,25 (6) 
i=l 


where C;,(Pg;) is the fuel cost of the ith generator, P;; is the active power generated by the thermal 
generators, a;, b; and c; are the cost coefficients of ith generator. 
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— Equality constraints: 
The equality constraints represent the flow equations of the balanced powers as follows: 


N 
P, -P, =V,>V,(g,, cosd, +z, sind, ) (7) 
jal 
N 
O, - 2, =ViDV;(g, sind, +z, cosd; ) (8) 
j=l 


— Inequality constraints: 
The equality constraints represent the limits of variable control and state control of the power system and can 
be given as follows: 


min max 
Foi Ss PGi Ss F Gi 


min max 
0;; Ss Oz; S 0c; 

min max 
Vai S Voi S Vai 

min max 
Tri s Turi S Tyr 


min max 
SVCi Ss Osvci = SVCi 


max 
| sii = STi 


(9) 


The vectors of control variables u, and uz are respectively the cases without and with the presence 
of SVC devices on the power system, and can be described as follows: 


b= | Pe | (10) 


y= [Pezes Pepi V open Vesios Ovvens Ones | (1) 


Where: Pe are the active powers generated, V, is the generator voltage and Qcoyc is the reactive power 
injected by the SVC device. 


4. THE ANT LION OPTIMIZATION (ALO) ALGORITHM 

The ant lion optimizer (ALO) is considered as the most recent nature-inspired proposed by [21]. 
The modeling of the ALO algorithm based on the hunting mechanism of antlions in nature. The main 
objective of the ALO algorithm is to solve any optimization problems of constrained engineering, it can get 
an optimal solution for minimizing the objective function by satisfying various constraints. In the ALO 
mechanism, it can be hunting the prey (ant) through five main steps as follow; random walk of ants, building 
traps, trapping in antlions traps, sliding prey toward antlion and final step are catching preys and rebuilding 
traps for a new step of hunting. 

The ALO method mimics the hunting behavior of ant lions, the expression mathematically of the 
random walks of ants to detect the location of food is describes as follow: 


X (t) = [0, cumsu (2r(t, ) — 1), cumsu (2r (t, ) — 1), ..., CUMSU (2r (¢, ) — 1)] (12) 


where X denotes the random walks of ants, cumsu is the cumulative sum, t is the step of random walk, n is 
the maximum iterations and r(t) show the stochastic function and given as follows: 


r(t) _ i if rand >0.5 


‘ (13) 
QO if rand <0.5 
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where rand represents a randomly number uniformly distributed in the range of [0,1]. 
The Details of different steps describe the relationship between predators and preys in the ALO 
method are explained as follow: 


4.1. Random walk of ants 
In every step of optimization in the ALO algorithm, ants move randomly inside the boundaries of 
the search space based on the (14), the random walks of ants are normalized by using the following: 


(x!~a)+(0— 
(b,-4;) 


X; = +¢; (14) 


where the a;; b; denotes the minimum and maximum of random walk respectively. cf and df are indicated 
the minimum and maximum of ith variables at tth iteration. 


4.2. Trapping in antlions traps 
The random walks of ants are influenced by antlions traps and are modeled as follows: 


C= Antlion’, Sey od (15) 
d; = Antlion’, + d' (16) 


4.3. Building traps 
In this work, the ALO algorithm is required to use a roulette wheel selection operator for selecting 
the better antlions based on their higher fitness, forgive a high chance for catching ants. 


4.4. Sliding ants toward antlion 

When the ants move toward near the center of the pit. However, once antlions realize that an ant is 
in the trap, they shoot sands outwards the center of the pit. To model this mechanism mathematically, the 
radius of the ant’s random walk is decreased correspondingly using (17) and (18): 


c aa (17) 
I 
d' 

ee 18 
F (18) 


4.5. Catching preys and rebuilding the traps 

The final step of hunting is when the prey reaches into the bottom of the antlion pit and is caught in 
the antlion’s jaw. After this stage, the antlion pulls the prey inside the sand and consumes its body. Then the 
antlion updates its new position to the latest position of the ant, to enhance its chance of catching new prey. 
The equation which models the catching prey and rebuilds the pits is given as follows: 


Antlion; = Ant), if f (Ant!) > f (Antlion' ) (19) 


where Antlion{, Ant! represents the position of the selected jth antlion and ith ant at iteration t. 


4.6. Elitism 

The elitism of an ant lion is determined by using the roulette wheel selection (RWs) at each step of 
optimization. The best antlion selected should be capable to affect the movements of all the ants 
at any iteration is saved as elite. The elitism mechanism for repositioning of a given ant described in the 
following equation: 


Ri +R 
Ant, = “ae (20) 
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where Rj§ is the random walk around the selected antlion using the roulette wheel at t — th iteration, Rf is the 
random walk around the elite antlion at t — th iteration. 


5. RESULTS AND ANALYSIS 

In order to show the performance and efficiency of the proposed algorithm ALO to solve 
optimization in the larger system dimensions, the OPF has been performed on the Algerian 114-bus system. 
This system network involves 15 generators, 175 lines, 16 tap changer transformers are located from line 160 
to line 175 and 99 load bus of total demand are 3,727 MW and 2070 MVar. The economic and technical 
parameters of 15 generators in the of the Algerian 114-bus power system in Ref. [22]. In this study, the 
proposed algorithm has been applied on the system under two different simulation cases that are considered, 
with and without static VAR compensator (SVC) devices. In the two simulation cases, 30 independent runs 
were executed for establishing the superiority of the ALO method with the population size equal to NP = 40 
and the maximal iterations are 200. The flowchart of the implementation of the proposed algorithm for 
solving the OPF problem that minimizes the total generation cost is shown in Figure 1. 


- ™ 
Start 


Read System Data Algerian 114-bus power Using (14) for create a random walk for the i'* ant 
svstem and normalize it 


Set the parameters of ALO method Update the position of ant by using (20) 


Initialize position of ant and antlions using — 


uniformly distributions , _-—~ Variables ~~ No, Tagging at 
Sr the limits 
"Yes 


Map control variables from ants into load flow 


Calculate fitness values of all ants 


Evaluate the TGC of the initial population 
from load flow 


Replace an antlion with its corresponding 


ant if it becomes fitter using (19) 


Update the elite if an Antlion becomes better 
Best results obtained = Best TGC 


Figure 1. Flowchart of the proposed algorithm applied to optimization 


No Iter < Itetmay 


Yes 


Select an antlion (using the roulette wheel) 


Slide ants towards the antlions 
using (17) and (18) 


5.1. OPF for Algerian electrical network system without SVC devices 

In this first case, we perform simulations on the Algerian 114-bus system without the SVC device. 
The objective function used is minimizing the fuel cost of 15 thermal generators and the vector of control 
variables contains the active powers generated as shown in (10). The optimization results obtained by the 
ALO algorithm compared with the grey wolf optimizer (GWO) and other optimization methods in the 
literature are tabulated in Table 1. The corresponding convergence of the proposed algorithm and GWO 
algorithm is shown in Figure 3. 
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Table 1. Comparison of TGC obtained by ALO algorithm and other algorithms in the literature 


Variables eae ALO Gwo DE [22] ra QP [24] | MOALO [25] 
Po(MW) 135 1350 438.2927 449.2373. 462.3008 ~~~ 455.9113 449.559 458.0600 
Pos(MW) 135 1350 452.1957 455.3286 459.5589 455.9219 449.559 451.1905 
P51,(MW) 10 100 99.8645 99.9064 99.9431 100.0000 100.000 74.91904 
Pois(MW) 30 300 199.3340 193.9445 «192.5196 «194.3179 195.368 212.0149 
Po,(MW) 135 1350 444.4851 453.0225 453.0142 «448.7254 449.559 436.8922 
Pe(MW) 34.5 345 204.8873 193.5834 196.6569 —«:196.0150 195.368 236.6123 
Poox(MW) 34.5 345 207.6710 191.9811 189.0239 190.8388 195.368 197.7529 
Posx(MW) 34.5 345 192.2261 185.6850 193.9372 _—:197.8609 195.368 246.6429 
Pooo(MW) 34.5 345 194.8395 193.3146 192.1215. —«*193.7858 195.368 171.7571 
Pras (MW) 30 300 191.3274 189.5693 «188.1283 ——-190.9545 195.368 163.4842 
Pog (MW) 30 300 184.6229 196.7435 «189.0847 «191.9255 195.368 214.0323 
Poroo(MW) 60 600 600.000 600.0000 599.9752 600.0000 600.000 599.9999 
Poror(MW) 20 200 200.000 200.0000 199.9703 200.0000 200.000 199.9999 
Peroo(MW) 10 100 99.9996 100.0000 99.9909 100.0000 100.000 67.19996 
Poii(MW) (10 100 93.5987 99.8715 99.9415 100.0000 100.000 81.96231 
Total generation cost ($/h) 19141.7714  19171.9582 9203.34 19199.444. _19197.696 _19355.859 
Active power loss (MW) 76.3446 75.1879 89.2570 89.2570 89.257 - 


The optimal values of control variables, total generation cost, and active power losses are 
summarized in Table 2. Over 30 independent trial runs were executed in this case as shown in Figure 2. 
From this results obtained, the best value of total generation cost and active power losses by the proposed 
method are 19141.7714 $/h and 76.3446 MW respectively, these values are better than the results obtained 
by different algorithms previously reported in Table 1. The Figure 3, allows us to note, in the first place, 
that the method ALO converges towards the global optimum at the iteration 100 while the convergence 
of the GWO method is reached at the iteration 180. So, the results obtained showed the proposed method 
ALO superior and robust compared to the GWO method in terms of getting the best solution for solving the 


OPF problem. 


19600 
19500 
19400 
19300 
19200 
19100 


Fuel cost ($/h) 


Trial run number 


Figure 2. Distribution of minimum total generation cost obtained for 30 trails by ALO 
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Figure 3. The Convergence of ALO and GWO for the Algerian 114-bus power system 
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5.2. OPF for Algerian electrical network system with SVC devices 

In this case, the SVC devices are implemented in the Algerian 114-bus system to improve the 
voltage profile and reduce the TFC. The vector of control variables, in this case, contains the active powers 
generated, the generator voltage and the reactive power injected by the SVC devices as shown in (11). 
The first optimal placement of the SVC device in the Algerian 114-bus system at bus N°89 (Souk Ahras), 
and bus N°68 (Sedjerara). The optimization results are given in Table 2. 
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Table 2. Optimization results obtained by ALO with different optimal placement of SVC 


Variables Limites _ Without SVC SVC Bus N°89_ SVC Bus N°89_ SVC Bus N°68 & 89 
Min Max 
Pg4a(MW) 135 1350 439.6459 428.7504 422.0728 432.5625 
Pos(MW) 135 1350 448.1696 435.1449 426.7332 439.0722 
Po, MW) 10 ~=100 82.0231 93.1780 100.0000 100.0000 
Pois(MW) 30 =300 171.8617 223.7201 235.9225 229.0151 
Poi7(MW) 135 1350 443.1784 437.5715 420.1050 441.9554 
Porg(MW) = 34.5345 209.4855 196.8190 184.4021 162.6831 
Po22(MW) 34.5 345 186.9316 187.4384 167.9495 218.4372 
Pos2(MW) 34.5 345 216.0484 224.8141 221.2963 198.0932 
Pogo (MW) 34.5 345 232.0432 211.1695 198.3792 180.6525 
Pog3(MW) 30 = 300 192.5640 185.8489 224.2547 221.0994 
Pcog(MW) 30 300 —:165.8376 174.1723 186.9096 161.7763 
Peroo(MW) 60 600 599.9998 600.0000 600.0000 600.0000 
Poi91 MW) 20 200 199.9999 200.0000 200.0000 199,9999 
Po199 MW) 10 100 99.9999 94.2300 99.2689 99.9999 
Poi11 (MW) 10 +100 99.9997 94.3467 99.9996 100.0000 
Vmin 1 114 0.9023 0.9022 0.9273 0.9393 
Qsvc68 (Mvar) -45 45 - 20.9076 - 20.9159 
Qsvc89 (Mvar) -45 45 - - 33.4490 33.3790 
Total generation cost ($/h) 19061.4915 19042.7382 19023.2724 18999.6809 
Active power loss (MW) 60.7884 60.2037 60.2935 58.3467 


Table 2 shows the results of TFC and real power losses obtained by using the proposed algorithm in 
the case without and with SVC devices, separately or multiple in buses N°68 and N°89. From these results 
obtained, it can be observed that the presence of SVC devices in all cases improved considerably the TGC 
and active power loss. The convergence curve of the proposed algorithm is shown in Figure 4. From this 
figure, we notice that the algorithm ALO converges towards the global optimum at the iteration 100 for all 
cases study when the SVC devices installed. 
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Figure 4. Convergence plot of ALO methods in the Algerian 114 bus power system 


On the other hand, SVC devices have managed to improve the voltage profile as shown in Figure 5, 
from this figure, it can be noted that the case without the presence of SVC, the two circles in this figure 
determine the two areas that have critical load bus voltage in Algeria 114-bus system. So, when the SVC 
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devices were installed at busses 68 and 89 separately or multiple, the voltage values of this optimal 
emplacement were increased to | p.u as shown in Figure 5, respectively. Furthermore, this increment in the 
voltage values in the optimal placement of SVC allows for improving the voltage at the critical load buses 
compared to the previous case (without the presence of SVC device). 


1,15 


Without SVC m SVC N°68 m SVC N°89 SVC N°68 & 89 


Voltage profile (pu) 
— 


0,9 


0,85 


Figure 5. The effect of SVC device on the voltage profile in the Algerian 114-bus power system 


6. CONCLUSION 

In this paper, we have validated the new metaheuristic technique, called, Ant Lion Optimizer (ALO) 
for real and large scale Algerian 114-bus power system to solve optimal power flow (ORF) problem. The 
ALO algorithm was successfully applied to solve the OPF problem with and without SVC devices. From the 
results obtained in the case without the SVC device, the proposed algorithm has been the best result 
compared with the method developed by us, called, grey wolf optimization and other methods in the 
literature defined in this paper, like DE, GA-ED-PS, and OP. In the case with the presence of SVC devices, 
the ALO algorithm was used to identify the optimal sizing and placement of SVC devices in the Algerian 
114-bus system based on the location of the lowest voltage load buses in the power system. The optimization 
results achieved by using the ALO algorithm with presence the of SVC devices given the best results to 
minimize the total fuel cost, reduce the active power losses and improving the voltage profile based on the 
optimal placement and sizing of SVC devices. Based on the results of both case studies in this paper, it can 
be concluded that the ALO algorithm is capable of solving the OPF problem for a large scale power system 
with and without the presence of SVC devices. 
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